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We present new results in an ongoing study of the nature of the high temperature crossover in QCD with two 
hght fermion flavors. These results are obtained with the conventional staggered fermion action at the smallest 
lattice spacing to date — approximately 0.1 fm. Of particular interest are a study of the temperature of the 
crossover, an important indicator of continuum scaling, a determination of the induced baryon charge and baryon 
susceptibility, used to study the dissolution of hadrons at the crossover, the scalar susceptibility, a signal for the 
appearance of soft modes, and the chiral order parameter, used to test models of critical behavior. 



1. INTRODUCTION 

Lattice simulations of high temperature QCD 
have provided our only present firmly grounded 
theoretical insights into the phenomenology of 
the crossover from hadronic matter to the quark- 
gluon plasma and into the nature of the plasma. 
Outstanding problems include (1) establishing 
continuum scaling, important not only for deter- 
mining the temperature of the crossover, but nec- 
essary for the validity of all dynamical fermion 
simulations, (2) establishing the mechanism for 
the dissolution of hadrons at the crossover, (3) 
exploring the intricate critical behavior associ- 
ated with the phase transition, and (4) obtain- 
ing a quantitative characterization of the quark- 
gluon plasma, including the equation of state. 
To achieve these goals requires a combination 
of advances in algorithms and computing power 
[1]. Rapid improvements in lattice formulations 
may soon have a significant impact on dynamical 



fermion simulations [2]. Here we give a progress 
report on an analysis of simulations with the con- 
ventional staggered fermion action at the smallest 
lattice spacing to date, namely, approximately 0.1 
fm [3]. 

The simulation with two fiavors of quarks was 
carried out at two quark masses, ma = 0.008 and 
0.016 and six couplings, = 5.65, 5.70, 5.725, 
5.75, 5.80, 5.85, except that the 5.725 coupling 
was not simulated at the higher mass. In each 
case the simulation was extended to at least 2000 
molecular dynamics time units. Lattices were 
saved at intervals of 8 time units. For the present 
analysis the first 500 time units were omitted. 
Most of the results reported here are based on 
an analysis of the approximately 180 remaining 
lattices at each parameter pair. Although results 
presented here are based on more data than re- 
ported at the conference, they should still be re- 
garded as preliminary. 



*Presented by C. DeTar 
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Figure 1. Chiral condensate {ipip) vs Q/g"^ 
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Figure 2. Polyakov loop i)s 6/(/^ 



2. LOCATING THE CROSSOVER 

2.1. Polyakov loop, {ipip), fuzzy loop 
Plotted in Figs. 1 and 2 are the observables 

traditionally used to locate the crossover, which 
is signalled by an inflection point. Despite small 
errors, it is evident that locating an inflection 
point from these results is challenging, to say 
the least. To strengthen the crossover signal we 
tried constructing a "fuzzy" loop variable by re- 
placing each time- like link Ut in the conventional 
Polyakov loop with 

t^fuzz.t = aUt + ^/staple (1) 

where [Staple is one of the six staples of length 
three links connecting the ends of the conven- 
tional single link [4]. This deflnition is not recur- 
sive and the resulting link is not projected onto 
SU(3). For present purposes we chose a weighting 
a = l3 = 1/7. Although this variable is slightly 
less noisy than the Polyakov loop, the crossover 
signal remained unclear. 

2.2. Baryon susceptibility 

The conventional Feynman path integral sim- 
ulates the grand canonical ensemble in baryon 
number at zero chemical potential. The baryon 
susceptibility measures fluctuations in the baryon 



number of the ensemble. It is deflned as the 
derivative of the baryon charge density with re- 
spect to chemical potential. The susceptibility 
can be deflned separately for each flavor. Thus 
with two quark flavors, two susceptibilities can 
be measured: a flavor singlet and flavor nons- 
inglet [5]. Both quantities are obtained at zero 
chemical potential. 

We present here only Xns- Results are com- 
pared with those for lower Ni in Fig. 3 [5]. Also 
indicated is the free lattice quark value for each 
Ni. A common pattern is the tendency for the 
susceptibility to rise abruptly at crossover and ap- 
proach the free lattice quark value asymptotically. 
At lower values of Ni, where the crossover has 
been located with traditional methods, we flnd 
that the baryon susceptibility reaches 1/3-1/2 of 
the free quark asymptotic value at the crossover. 
Following the same rule for Ni = 12 places 
the crossover in the ma = 0.008 series between 
= 5.65 and 5.70 and in the ma = 0.016 se- 
ries between G/g"^ = 5.75 and 5.80. 

To convert these results to a temperature, we 
use a scale in which the rho meson mass is taken 
to be 770 MeV, regardless of quark mass. The 
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Figure 3. Nonsinglet baryon susceptibility v.s6/g^ 
for Ni = 4,6,8, 12. Free lattice quark values for 
each Ni are indicated by horizontal lines on the 
right of the plot. 

rho mass in lattice units is obtained in turn from 
fitting a polynomial to masses in a compilation 
of two-fiavor staggered fermion spectral simula- 
tions [6], and includes an extrapolation beyond 
the parameter range of spectral measurements 
(Q/g^ > 5.7) using tadpole-improved asymptotic 
scaling [3]. From the baryon susceptibility we 
then place the crossover at Tg = 143 — 154 MeVp 
at the lighter quark mass and Tg = 142 — 150 
MeVp at the heavier quark mass. These values 
are consistent with results at smaller Ni [3], so 
reveal no drastic departures from scaling. 

2.3. Induced quark number 

Another confinement-sensitive observable is the 
induced quark number [7]. This observable mea- 
sures the total residual light-quark number in an 
ensemble containing a single test quark. 



Qmd = / Pind(r)d^r 



(3) 



The induced quark number density is measured 
in the presence of the test quark. Operationally, 
the quantity measures the correlation between the 
Polyakov loop and the light-quark density [7]. To 
reduce the considerable noise, we introduced the 
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Figure 4. Induced quark number. Octagons are 
for fixed lattice spacing. 

test charge through the fuzzy Polyakov loop vari- 
able described above. Results are shown in Fig. 4 
and compared with results from simulations at 
lower Ni. The induced quark number is expected 
to be exactly —1 at zero temperature, since con- 
finement requires screening of the test charge by a 
single antiquark. At the crossover, this quantity 
rises rapidly, approaching zero in the high tem- 
perature phase. At lower Ni the induced quark 
number reaches approximately —0.1 at crossover. 
Applying this rule to the Ni = 12 data gives 
= 5.65 - 5.70 (T = 143 - 154 MeV^) in 
the ma = 0.008 series and 6/g^ = 5.70 - 5.80 
(T = 134 - 150 MeYp) in the ma = 0.016 se- 
ries. This crossover location is consistent, but 
somewhat less precise, than that found from the 
baryon susceptibility. 

3. CHIRAL SUSCEPTIBILITY 

Another signal for the crossover is the chiral 
susceptibility 



d {ipip) 
dm 



(4) 



6/3 = 



which measures fiuctuations in the chiral order 
parameter (i/'i/') [8]. Since it is the space-time 
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Figure 5. Disconnected chiral susceptibility for 
Nt = 4[8] and 12. 

integral of the correlator (^tptp(r)tptp(O)'^ , a peak 
in this observable occurs at a minimum in the a 
meson screening mass, indicating the presence of 
a soft mode. Such soft modes are expected in 
models of critical behavior [9]. Results for the 
quark-line disconnected contribution to this sus- 
ceptibility are compared in Fig. 5 with the results 
of Karsch and Laermann at = 4 [8]. The trend 
is consistent with a peak at the = 5.70 data 
point (153 MeVp) in the lighter mass data and at 
the6/g^ = 5.80 data point (150 MeVp). However, 
the signal in the higher mass data is not strong 
enough to be predictive. That the signal should 
be weaker at higher quark mass is expected from 
models of critical behavior that place the critical 
point at zero quark mass. The mass ratio m^/m^ 
measures proximity to the critical point. In the 
Ni = 4 data, this ratio is less than 0.1 for the 
ma = 0.020 data and between 0.1 and 0.2 for the 
ma = 0.0375 data. For the Ni = 12 data, this 
ratio is between 0.2 and 0.3 for ma = 0.008 and 
0.3 to 0.4 for ma = 0.016. 

4. CONCLUSIONS 

With the conventional staggered fermion ac- 
tion, locating the high temperature crossover at 
0.1 fm for m^/m^ > 0.2 is difficult in the con- 
ventional observables (Re_P) and {ipip) . How- 



ever, the baryon susceptibility and induced quark 
number provide an adequate determination. The 
crossover temperature at the lighter quark mass 
is found to be in the range 143 — 154 MeVp, con- 
sistent with results at lower Nf. Other qualita- 
tive features of the crossover, seen at lower Ni, 
are confirmed. In particular, the crossover is sig- 
naled by an abrupt rise in baryon susceptibility 
and induced quark number. 
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